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We study theoretically many-body equilibrium magnetic phases and corresponding thermody¬ 
namic characteristics of ultracold three-component fermionic mixtures in optical lattices described 
by the SU(3)-symmetric single-band Hubbard model. Our analysis is based on the generalization 
of the exact diagonalization solver for multicomponent mixtures that is used in the framework of 
the dynamical mean-field theory. It allows us to obtain a finite-temperature phase diagram with 
the corresponding transition lines to magnetically ordered phases at filling one particle per site (1/3 
band filling) in simple cubic lattice geometry. Based on the developed theoretical approach, we 
also attain the necessary accuracy to study the entropy dependence in the vicinity of magnetically 
ordered phases that allows us to make important predictions for ongoing and future experiments 
aiming to approach and study long-range-order phases in ultracold atomic mixtures. 
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I. INTRODUCTION 

There has been remarkable progress in approaching 
quantum magnetism in ultracold two-component atomic 
mixtures in optical lattices, where short-range antifer¬ 
romagnetic correlations were detected and measured re¬ 
cently Qi- It is believed that observations and sub¬ 
sequent detailed studies of the related long-range-order 
states of matter, together with other novel quantum 
phases, will help not only in designing new generations of 
strongly correlated materials, but also in quantum sim¬ 
ulation applications. At the same time, the mentioned 
experimental success, together with still existing chal¬ 
lenges (e.g., limitations in cooling Q), shifts the research 
horizons further. 

From a theoretical point of view, while the main 
ingredients for long-range-order magnetic phases (e.g., 
phase diagrams and critical entropies) are determined 
with a high accuracy for the single-band two-component 
fermionic Hubbard model (see Refs. ii for details), 
there are many fewer quantitative predictions for mul¬ 
ticomponent mixtures. However, there is clear evidence 
that these systems possess a very rich physics that re¬ 
sults in a variety of quantum many-body phases and their 
unique properties |6l4l2l|. 

In this paper we present a theoretical approach that 
allows us to study long-range-order magnetic phases in 
the SU(3)-symmetric Hubbard model in more detail and 
determine the finite-temperature phase diagram struc¬ 
ture depending on the main system parameters, the hop¬ 
ping amplitude and the interaction strength, that can 
be tuned with a high degree of freedom for ultracold 
atoms in optical lattices. However, we not only aim to 
bridge the mentioned gap, but also analyze by this exam¬ 
ple whether multicomponent mixtures can be considered 
advantageous in comparison with their two-component 
counterparts for the purpose of approach ing spe cific long- 
range-order states as suggested in Refs. [la - H^ . 


II. SYSTEM AND MODEL 

We consider a system described by the SU(3)- 
symmetric Hubbard Hamiltonian of the type 
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where t is the hopping amplitude between nearest- 
neighbor sites i and j with the corresponding shorthand 
notation {ij) under the sum, U is the on-site interac¬ 
tion strength between atoms in different internal (e.g., 
hyperfine) states a and /3, {cia) is the fermionic cre¬ 
ation (annihilation) operator of the atom in the state a 
on the lattice site i, and hia = cJ^Cia and hi = 
denote operators of the particle number. In the second 
line, we also introduced the amplitude of the external 
trapping potential Vi on the lattice site i and the chem¬ 
ical potential /i that sets the total number of atoms in 
the system. We assume that the optical lattice is deep 
enough Idat ^ 5Er {Er is the recoil energy of atoms) for 
the single-band approximation to be correct. 

The system described by the Hamiltonian m can be 
experimentally realized by loading ultracold atoms that 
are prepared in three different quantum hyperfine states 
to the optical lattice. Since atoms in all three internal 
states (called spin components below) must interact with 
each other with the same strength U, appropriate can¬ 
didates are mixtures of alkaline-earth atoms (^^^Yb and 
®^Sr; see, e.g., Refs. Bill) and mixtures of atoms, 
where by means of a proper choice of hyperfine states and 
Feshbach resonances one can realize a regime with almost 
equal scattering lengths between spin components. In 
both experimental cases the lattice depth can be easily 
tuned by the intensity of laser beams, thus the coupling 
Uft can be changed in a wide range. 
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It is worth emphasizing that in the case of a homoge¬ 
neous system (V^ = 0) and filling one particle per site 
(i.e., n = 1 or, equivalently, 1/3 band filling that is, 
strictly speaking, only approximately fulfilled by setting 
^ Ri {7/2; see Ref. [l^ for details) in the strong-coupling 
limit U/t:^ 1 and by using an established technique [l^ 
one arrives at the SU(3)-symmetric Heisenberg model 
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with positive (antiferromagnetic) exchange coupling J = 
At'^/U and the local pseudospin projection operators 
Ski = ^cl^i^kapcpi, which are expressed in terms of 3 x 3 
unitary Gell-Mann matrices Xk = {Ai,..., Ag} [l^. Be¬ 
low, we do not limit ourselves to the strong-coupling 
limit; we study thermodynamic properties of the Hub¬ 
bard model (H]) at filling n = 1 depending on the coupling 
U/t and the temperature T/t. At the same time, the ef¬ 
fective Hamiltonian ([2]) is useful for a proper analysis of 
particular effects obtained by the theoretical approach 
introduced below. 


III. METHOD 

We use a dynamical mean-field theory (DMFT), the 
general idea and specific details of which can be found 
in Ref. [2^. Regarding other state-of-the-art numeri¬ 
cal methods, it is important to note that DMFT can 
be considered at present as one of the most optimal ap¬ 
proaches for the system under study. In particular, the 
determinant quantum Monte Carlo method [2l| cannot 
be applied for the SU(3)-symmetric mixture due to the 
presence of the sign problem that cannot be suppressed 
or excluded as, for example, in the SU(4)-symmetric mix¬ 
ture at half filling [l^. At the same time, the direct exact 
diagonalization approach is very limited in the number 
of lattice sites for such systems @ that may cause sig¬ 
nificant difficulties in a reliable finite-size scaling analy¬ 
sis for magnetic phase diagrams. Therefore, below, we 
briefly introduce key generalizations and main steps in 
DMFT that allow us to obtain the central results of the 
paper and can be important for further studies of mul¬ 
ticomponent mixtures in the framework of the approach 
used. 


A. Generalization of the exact diagonalization 
impurity solver for multicomponent mixtures 

The central idea of DMFT consists in mapping the 
original lattice problem © that is usually intractable to 
the local (impurity) problem that can be solved exactly. 
Taking the Anderson impurity model (AIM), one can no¬ 
tice that it can be extended to mixtures with a multiple 
number M of spin components in a straightforward way. 
Since we use the exact diagonalization (ED) solver in our 


simulations (originally introduced in Ref. [2^ for two spin 
components A/” = 2), the corresponding impurity Hamil¬ 
tonian can be written as 
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where a, a' = {!,..., A/”} denote the spin indices and 
the index I = {2,..., Us} labels the number of the bath 
orbitals in the AIM with rig being the cutoff number pe¬ 
culiar to the ED approach. The two terms in the first 
line of Eq. Q correspond to the energies of electrons in 
the bath and the hybridization between the bath and the 
impurity, respectively. The operators {aia) and dj. 
(do-) are the creation (annihilation) operators of electrons 
on the bath’s orbital I and the impurity, respectively; 
the quantities e/o- and Viu are the so-called Anderson pa¬ 
rameters that set the amplitude of the processes in this 
model. These parameters are determined iteratively in 
the DMFT approach. The terms in the second line have 
a direct correspondence to the original Hubbard Hamil¬ 
tonian m as a generalization to the SU(A/')-symmetric 
mixture. 

In the Fock-space representation, the occupation of the 
orbital I by the spin component a can take two values 
nia = {0,1}, thus, in the most general case, the matrix 
that is necessary to diagonalize within the ED solver is 
of the size Lq x Lq with Lq = 2"'’^. This results in a 
strong limitation of the number of orbitals rig that is fea¬ 
sible to account in numerical calculations for the given 
number of spin components Af. However, analogously to 
Refs. [ 20 , 123 , due to the structure of the Hamiltonian Q 
that corresponds to the original lattice problem o , we 
note that in the case under study it does not mix differ¬ 
ent spin sectors, |nii,..., • • • \nij\f ..., n„„ 7 V'), thus 

the total charge per spin component q,j = '^i <7 can 

be considered as a conserved quantum number. Defin¬ 
ing the quantum state as \Qi) = |(;|*\ ..., we de¬ 
termine in that way the total number of distinct states, 
z = {!,..., (ng -h 1)^}. 

Therefore, the matrix Ap x Ap diagonalization problem 
transforms to diagonalizing (rig -1-1)^ blocks of a differ¬ 
ent (but much smaller) linear size Li expressed in terms 
of the binomial coefficients Li — H^i "^he de¬ 

scribed procedure allows us to account for more orbitals 
in the impurity problem, thus increasing the correspond¬ 
ing accuracy of results. In particular, from the direct 
calculation analysis, we conclude that by comparing to 
two-component mixtures with rig = 7, it is completely 
feasible to use rig = 5 for three-component and rig = 4 for 
four-component mixtures within the ED solver in DMFT 
approach. 

After the diagonalization is performed, we obtain the 
interacting Green’s functions in the Matsubara-frequency 
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space according to the standard definition (see 


also [ 2 ^) 
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where Z = e is the partition function, the sum 

is taken over the full set of eigenstates \j) and \k) with the 
corresponding energies Ej and Ek^ujn = 7r(2n+l)T is the 
fermionic Matsubara frequency, and T is the temperature 
(we use units such that ks = !)• Subsequently, the local 
self-energies for each spin component are calculated from 
the Dyson equation 

Y.a{i0Jn) = Qa^{i0Jn) - G~^{iUJn), ( 5 ) 


where the Weiss Green’s function is defined by the set 
of Anderson parameters entering the Hamiltonian (|5|) 
through the analytic relation Q~^(iu}n) = + ^ — 


B. Clustering in the real-space lattice projection 
for homogeneous systems 


To obtain the lattice Green’s functions for self- 
consistency conditions in DMFT, one can use the non¬ 
interacting density of states D{e) for a particular lattice 
geometry and the local self-energies obtained within 
the impurity solver. However, this method can be suc¬ 
cessfully applied (see, e.g.. Ref. for more details when 
there is no translational invariance broken (e.g., in the 
paramagnetic or ferromagnetic phase). 
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or the symmetry is broken towards a bipartite structure 
(e.g. charge-density wave or antiferromagnetic Neel-type 
ordering) 


Gas 
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where (as = + k- ~ ^asi'ii^n), the sublattice indices 

s = A,B, and its opposite s = B,A. 

Obviously, there is no general expression possible in 
terms of D(e) for more exotic types of magnetic or¬ 
der and, in particular, for the three-sublattice antiferro¬ 
magnetic state that is peculiar to three-component mix¬ 
tures at low temperatures [1, [l^. Therefore, a real- 
space generalization of DMFT (RDMFT) (introduced in 
Refs. [^) can be considered as a proper extension 
to account for more types of magnetic order in an unbi¬ 
ased way [13 • Within this approach, the corresponding 
lattice Green’s function is obtained from the inversion of 
the following real-space matrix: 

(^^n)]ss^ — E k (^a^n)]^ss' GsS (^) 


where the indices s, s' denote the lattice sites, 6ss' is the 
Kronecker symbol, and is the hopping matrix element 



FIG. 1. (Color online) Schematic illustration of the clustering 
procedure in RDMFT by an example of the three-sublattice 
structure in a square lattice geometry {L is the number of 
sites in the x direction) that allows for the antiferromagnetic 
phase with the ordering wave vector Q = (27r/3, 27r/3). 


that equals t if there is tunneling possible between the lat¬ 
tice sites s and s' and zero otherwise, which depends on 
the lattice geometry, system size, and the boundary con¬ 
ditions (open or periodic) used in the original model ([T|) . 

It is important to note that RDMFT requires more 
computational resources in comparison to the single-site 
® or two-sublattice ® DMFT. It is caused not only 
by the necessity of the matrix inversion (which is feasi¬ 
ble for systems consisting of sites or even more), 

but also because the self-energies T,as{iuJn) must be cal¬ 
culated within the impurity solver on every lattice site s 
of the original lattice problem O- However, the lat¬ 
ter limitation can be significantly reduced for homoge¬ 
neous systems in the case when the type of breaking of 
the translational symmetry is known beforehand (e.g., 
from unbiased RDMFT calculations for smaller system 
sizes). In Fig. [T] we show one particular type of clus¬ 
tering that allows for the possibility of three-sublattice 
antiferromagnetic ordering in the system. This simplifi¬ 
cation requires the impurity problem to be solved only 
on three lattice sites, while the Green’s functions are ob¬ 
tained from Eq. ([5]) for the full-size system, thus effec¬ 
tively minimizing the influence of the finite-size effects. 


C. Entropy analysis 

Concerning experiments with ultracold atomic mix¬ 
tures in optical lattices, a crucial quantity for approach¬ 
ing magnetically ordered phases is not the temperature 
(which is a natural variable in the DMFT approach op¬ 
erating in the framework of the grand-canonical descrip¬ 
tion) but the entropy, since atoms in optical lattices can 
be considered to some extent as isolated from the sur¬ 
rounding environment and key system parameters (e.g., 
the coupling strength U/t) can be changed adiabatically. 
Therefore, quantitative theoretical predictions for the en¬ 
tropy become of a high importance in this field of re¬ 
search. 

To calculate the entropy in our theoretical approach, 
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according to Ref. [2^, it is efficient to use the Maxwell 
relation for the entropy per lattice site ds/d^ = dn/dT. 
It can also be written in the form 


s(/r,C7,r) 



dn(fj,', U, T) 




(9) 


which is more suitable for the DMFT analysis, since /i, 
U, and T are input parameters, whereas the filling n is 
the local observable measured within the impurity solver. 

However, in the direct numerical analysis it is rather 
convenient to parametrize the chemical potential to the 
form fi{r) = ni — Vr^ that is analogous to the local- 
density approximation for our system in the external 
trapping potential of the parabolic shape with the am¬ 
plitude V and the radial distance r measured in units 
of the lattice constant. Therefore, now we also access 
the entropy distribution in the harmonic trap, since the 
entropy at the given radial point ri is defined as follows: 

dn 

-^rdr, ( 10 ) 

with the boundary condition at the edge point i?max of 
the trap s(i?i„ax) = 0 that is fulfilled in our calculations 
by the condition n(r > i?niax) = 0. 

It is worth noticing that for the entropy analysis of ho¬ 
mogeneous systems that is applied to the phase diagram 
with the fixed filling ni (e.g., ni = 1), the parameter 
ri in Eq. must be determined independently at ev¬ 
ery [U, T ) po int from the condition n(ri, U, T) = ni (see 
also Ref. [l2)- Probably, the only exception can be made 
in this respect for the case of half filling in the S\J{Af)- 
symmetric Hubbard model, where one obtains the fixed 
condition for the chemical potential /rhf = (A/” — l)t//2 
that is independent of quantum and thermal fluctuations 
due to particle-hole symmetry at any J\f. 


IV. RESULTS 

With the theoretical approach developed we study the 
magnetic (T — t/)-phase diagram. Concerning the spe¬ 
cific choice of the ED solver for this purpose in DMFT, 
we find that within the necessary accuracy it allows us 
to access low enough temperatures and has flexibility 
in the choice of the coupling strength UJt in the Hub¬ 
bard Hamiltonian m even for rig = 4. Here, according 
also to the limitations of the chosen impurity model ([3]) 
and analogously to Refs, [ill [13, we restrict ourselves to 
the case when the SU(3) symmetry can be spontaneously 
broken in two directions (i.e., along two easy axes that 
correspond to the natural-color basis) in the pseudospin 
space. In particular, we measure relative occupations by 
each spin component or, equivalently, analyze two local 
magnetizations on the lattice site s, 

= (Sas) = (nis - n2s)/2, 

m[^'> = (Sss) = (nis + n2s - 2n3s)/Vi, ( 11 ) 



FIG. 2. (Color online) Finite-temperature phase diagram for 
magnetic phases in the SU(3)-symmetric Hubbard model at 
filling one particle per site (n = 1) obtained by DMFT. At 
weak coupling and low temperatures (T < O.lt) the system 
size was chosen large enough (> 9® lattice sites) to minimize 
the influence of finite-size effects and to guarantee reliable 
numerical convergence of RDMFT. The reference data for the 
points Tci ,2 (CT-HYB) is taken from Ref. [l^ . 


and their periodic behavior in real space. Therefore, 
magnetically-ordered states are determined from the con¬ 
vergence analysis for the Green’s functions with differ¬ 
ent self-consistency conditions (ED-®, which allow for 
different types of the sublattice structure (i.e., different 
antiferromagnetic phases). Hence, the main results are 
summarized in Fig. ® 

By using the same scale in Fig. [3] we also show the 
transition line obtained by DMFT for the two-component 
SU(2)-symmetric mixture. This comparison directly 
shows a complete suppression of the antiferromagnetic 
order at weak coupling that we attribute to the absence of 
the key mechanism (leading to the so-called Slater-type 
antiferromagnet; see, e.g.. Ref. M) caused by the ab¬ 
sence of the perfect nesting of the Brillouin zone, since the 
band is only 1 /3 filled. In the opposite case of strong cou¬ 
pling U/t^ 1, as we mentioned, there is a direct map¬ 
ping to the effective spin model (®, thus the Heisenberg- 
type antiferromagnetic states arise at finite temperature 
in a cubic lattice geometry. This mechanism is also con¬ 
firmed in the phase diagram by a corresponding decrease 
of the critical temperatures rci, 2 , which are proportional 
in both cases to the magnetic coupling J = /U in the 
limit U/t ^ 1. 

Concerning the specific spatial structure of the ordered 
states, we show their possible in-plane configurations in 
Fig. 3 (for the sake of simplicity we omit the third spatial 
dimension in illustrations since the corresponding exten¬ 
sion is rather straightforward; see also Ref. [13 )■ Note 
that other equivalent spatial arrangements for possible 
three-sublattice antiferromagnetic (3s-AFM) states can 
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FIG. 3. (Color online) Schematic illustration of the observed 
magnetic phases at a filling of one particle per site. The color 
background indicates the presence of the additional net mag¬ 
netization in the 2s-AFM states. 


be obtained by 7r/2 rotations and spatial translations by 
the lattice constant along main axes. As for the two- 
sublattice antiferromagnetic (2s-AFM) states, while they 
allow for two possibilities in spatial arrangements in sim¬ 
ple lattice geometries, there is greater freedom for the 
symmetry breaking in the pseudospin space. In particu¬ 
lar, as shown in Refs. [Ill, fl^. for the SU(3)-symmetric 
mixture we observe two types of antiferromagnetic or¬ 
dering: a color-density wave (CDW2) and color-selective 
antiferromagnetic (CSAF) state. In particular, the spon¬ 
taneous breaking of the symmetry in the direction of 
the third spin component (again, restricting in measure¬ 
ments to easy axes) means that the CDW2 and CSAF 
states are characterized by = 0 with ^ 0 

and rn'^ = —m'g ^ 0 with = m'g > 0, respec¬ 
tively; see also Eq. (fTTI) and Fig. H] 

It should be mentioned that in our analysis of the 2s- 
AFM states we also observe the presence of the net mag¬ 
netization along the direction chosen by the symmetry¬ 
breaking mechanism. Moreover, from Fig.|4]one can con¬ 
clude that it has opposite signs in the CDW2 and CSAF 
states (i.e., n^A+n^B > 2/3 and n^A+n^B < 2/3, respec¬ 
tively) and increases with the temperature decrease. This 
means that at low temperatures the SU(3)-symmetric 
mixtures with the fixed equal number of atoms in three 
different hyperfine states can become unstable towards 
phase separation into the CDW2- and CSAF-ordered do¬ 
mains that compensate for the residual net magnetization 
produced by each other. Note that, according to Fig. SI 
the boundary between 2s-AFM and paramagnetic (PM) 
states in Fig. [2] (the line Tc 2 ) corresponds to the second- 
order phase transition similarly to two-component mix¬ 
tures. 

Concerning other phase boundaries in Fig. SJ it is nec¬ 
essary to note that the transition lines are obtained from 
the corresponding analysis of the magnetizations (ITT]) . At 
the critical temperature Td we observe discontinuities in 
these observables, i.e., the first-order phase transitions 
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FIG. 4. (Color online) Temperature dependence of the rela¬ 
tive occupation of two sublattices by spin components for two 
possible states in the 2s-AFM phase: CDW2 (left) and CSAF 
(right) at the fixed interaction strength U/t — 12. 


from the 3s-AFM state to both PM and 2s-AFM states. 
It should be mentioned that, regarding the explicit de¬ 
termination of the phase transition line in this case, it is 
necessary to analyze the grand-canonical potentials and 
thus study the system in more detail for possible coex¬ 
istence regions and metastable solutions. Here we must 
stress that in the case of a nonbipartite sublattice struc¬ 
ture, this task becomes highly nontrivial, thus transform¬ 
ing into a separate problem. However, the corresponding 
analysis of the impurity problem ([3]) can be performed 
within the ED solver with no significant effort since one 
has direct access to all eigenstates and partition func¬ 
tions that are also used in Eq. O- Hence, from these 
estimates we can conclude that there is no evidence of 
the coexistence between PM and 3s-AFM phases, but, at 
the same time, there are signatures that a narrow coex¬ 
istence region (with a height of AT ~ 0.02t) between 2s- 
and 3s-AFM phases can be present. Therefore, the actual 
transition line for Td can be expected to be lower from 
the side of the 2s-AFM phase than it is depicted in Fig. [2] 
(i.e., the line shown is effectively the upper boundary for 
the 3s-AFM phase). Let us emphasize that for a correct 
determination of the transition line one must solve a sep¬ 
arate problem and analyze the grand-canonical potentials 
corresponding to the original lattice Hamiltonian m in 
different phases with a proper account of the sublattice 
structure. 

Finally, we study the entropy of the system depending 
on the coupling strength Ujt and the temperature T jt 
at n = 1. This allows us not only to estimate the critical 
entropies necessary for approaching the magnetically or¬ 
dered phases specified in Figs. [2] and [3l but also to study 
an increase in the magnitude of the Pomeranchuk effect 
peculiar for the multiflavor mixtures. For this purpose, in 
Fig. [S] we also show the isentropic lines and the phase di- 
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FIG. 5. (Color online) Isentropic lines and magnetically or¬ 
dered phases of the SU(2)- and SU(3)-symmetric Hubbard 
models obtained by DMFT. 

agram for the SU(2)-symmetric mixture (the correspond¬ 
ing data are taken from Ref. [l^) that allow us to make 
a direct comparison between two- and three-component 
ultracold atomic mixtures in optical lattices from that 
perspective. 

From Fig. [5] we conclude that, despite the stronger 
suppression of the antiferromagnetic phases in units of 
the temperature Tjt^ the situation looks very optimistic 
for three-component mixtures from the point of view of 
the critical entropies necessary for approaching quantum 
magnetism based on the superexchange mechanism. In 
particular, there is a significant increase (approximately 
by 50%) of the critical entropy necessary for approaching 
2s-AFM states that we attribute to the stronger Pomer- 
anchuk effect. Of course, one should remember that 
DMFT provides inaccurate results for the values of the 
critical entropy in two-component mixtures due to its 
limitation on local correlations (see Ref. @ for more ac¬ 
curate data from exact methods). Naturally, the inac¬ 
curacy of the same origin should be expected in our es¬ 
timates for the SU(3)-symmetric mixture. However, the 
magnitude of the effect allows us to conclude that the 
advantageous properties should be confirmed with other 
theoretical methods and observed in direct experiments. 
According to our analysis, it is also a good sign that the 
critical entropy for the 3s-AFM phase < 0.5 is 

of a magnitude that is realistic to approach according to 
the current and expected progress in cooling techniques 
for ultracold fermionic mixtures. 


V. CONCLUSION AND OUTLOOK 

We developed a theoretical approach that allows us 
to study in detail the magnetic phase diagram of three- 
component fermionic mixtures in optical lattices with 
simple cubic geometry at finite temperature and a filling 
of one particle per site. At weak coupling for the SU(3)- 
symmetric mixture we observed a complete suppression 
of magnetic phases that arise then only in the region 


of the corresponding crossover to the Mott-insulating 
state. Within the two-sublattice and real-space DMFT 
generalizations used we studied two- and three-sublattice 
antiferromagnetically ordered states and analyzed their 
structure and critical temperatures. Despite the more 
complex structure of magnetic phases, in the entropy 
analysis we identified possible advantages of these mix¬ 
tures in comparison with the two-component counter¬ 
parts for approaching quantum magnetism in optical lat¬ 
tices due to more pronounced cooling that is facilitated 
by the stronger Pomeranchuk effect. 

From a theoretical point of view, many interesting 
questions are important to study further in this field of re¬ 
search. In particular, due to limitations of our approach, 
we were unable to make a strong statement regarding 
the quantum critical point at T = 0, i.e., whether it 
takes place at Uc ~ 9.6f (in accordance with the slope 
of the critical line in Figs. [2] and [S]) or there is an ex¬ 
ponential suppression of the critical temperature simi¬ 
lar to two-component mixtures. It seems that for this 
purpose DMFT with the corresponding generalization of 
the numerical renormalization-group solver to three- 
component mixtures that operates at T = 0 can be con¬ 
sidered as a good option. Another interesting direction 
for theoretical investigations is a more detailed analysis 
of coexistence regions and phase-separation effects both 
between 3s- and 2s-AFM states and inside the 2s-AFM 
phase (i.e., between CDW2 and CSAF states). 

Note also that in our calculations we restricted our¬ 
selves to the easy-axis directions in the eight-dimensional 
pseudospin space. Naturally, these configurations are not 
always preferred by the system consisting of ultracold 
atoms in different hyperfine states (see, e.g., Ref. for 
two-component mixtures). In particular, as mentioned 
in Ref. [^, in the case of breaking of the SU(3) sym¬ 
metry by population imbalance, additional terms appear 
(analogous to external magnetic fields along easy axes) 
in the effective Hamiltonian that push the system to¬ 
wards canted magnetic configurations. Therefore, this ef¬ 
fect must be properly accounted for in experiments (e.g., 
by introducing additional rf-pulse rotations in measure¬ 
ments of the magnetic phases studied). Note also that 
the easy-axis directions should be more feasible in sys¬ 
tems with the imbalance in hopping amplitudes. The 
latter can be successfully realized not only in mixtures of 
alkaline-earth atoms (^’’^Yb or ®^Sr) consisting of atoms 
in the ground and long-living metastable excited states, 
but also in mixtures, where a significant hopping 
imbalance can be effectively induced by lattice shaking 
with low heating rate [s^. Therefore, this would open 
a new direction in theoretical studies of possible advan¬ 
tages of the mass-imbalanced mixtures for the purpose of 
observations of magnetic phases studied in this paper. 

Another important direction concerns the inhomogene¬ 
ity and finite-size effects in these systems. It is stimulated 
not only by existing requirements from the experimental 
side, but also by theoretically studied effects of exotic 
Mott states 0' and spin separation of three-component 
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mixtures in the trap [l^ that could lead to finding new 
mechanisms in cooling of fermionic mixtures in optical 
lattices. Note that for the magnetically ordered states 
under study we assumed that the condition of one par¬ 
ticle per lattice site must be fulfilled. However, analo¬ 
gously to two-component mixtures, this condition does 
not need to be fulfilled exactly, since antiferromagnetic 
states are also observed at nonzero doping (e.g., the 3s- 
AFM state is stable for n £ [0.98,1.02] at U = 14t and 
T = 0.12f). Also, due to the fact that all the AFM states 
studied are accompanied by the Mott gap, there is a sig¬ 
nificant stability region for chemical potentials (e.g., the 
3s-AFM state is stable for ^ G [0.317,0.617] at 17 = lAt 
and T = 0.121:). The latter allows us to expect magneti¬ 
cally ordered domains of a sufficient size (with the linear 
size in one spatial direction of the order of ten sites or 
more, depending on the trap curvature) in the trapped 
systems, both theoretically and experimentally. 

Finally, the approach developed has the necessary re¬ 
quirements in accuracy for a detailed analysis of four- 
component mixtures. These systems are interesting, in 
particular, due to ongoing debates concerning theoret¬ 
ical determination of magnetically ordered states at a 


filling of one particle per site in the Hubbard SU(4)- 
symmetric model and the corresponding Heisenberg limit 
(see Refs. EinHii). Moreover, the detailed structure 
of the phase diagram at half filling is not completely de¬ 
termined [13 , thus the developed DMFT approach could 
be useful there as well. By means of four-component ul¬ 
tracold atomic mixtures in optical lattices, possibilities 
also arise to study from new perspectives the Kondo- 
lattice model as well as other multiband models with an 
analog of Hund’s coupling in solid-state materials. The 
latter can be realized now by tuning the relative ampli¬ 
tude of spin-exchange processes between particular spin 
components in these mixtures 
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